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ABSTRACT 

We present an analysis of the Fe II, Ge II, Mg I, Mg II, S I, S II, Si II and 
Zn II interstellar lines for 63 stars. The column density of Mg II, S II, Si II and 
Zn II is well correlated with the distance. However, the column density - distance 
relation should be used with care for the estimation of the distance to OB stars. 
For stars with large f{H2) this relation can lead to a large overestimation of the 
distance. Hydrogen - normalised column densities of Mg II, Si II and Ge II (our 
largest samples of data) drop with the interstellar reddening E(B-V) as expected 
for elements that are incorporated into dust grains. The Ge II abundance (Gell/H) 
is lower in dense molecular clouds. The abundances of all analysed elements are 
generally lower than their Solar System values. 
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1 Introduction 

McKee & Ostriker (1977) in their model have divided the interstellar medium 
(ISM) into three components: a low density hot ionised medium (HIM, n~ 10~^^ 
cm-^ T~ lO^ '^ K), a dense cold neutral medium (CNM, n~ lO^-^ cm'^, T~ 10^ '' 
K) and a photoionized cloud corona (T~ 8000 K). The cloud corona itself was 
divided into two subregions: warm ionized medium (WIM, fractional ionization 
xw/M ~ 0.7), and a warm neutral medium (WNM, ^wnm ~ 0.1). The filling factor 
of the cloud corona was estimated to be f~ 0.2, while the filling factor for the HIM 
was estimated to be f~ 0.7 0.8. The cold, dense clouds fill only f~ 0.02 0.04 of 
the interstellar space. Recently Cox (2005) wrote a review paper about the three- 
phase interstellar medium. 

Megier et al. (2005) pointed out, that the equivalent widths of the Ca II inter- 
stellar absorption lines are very well correlated with the distance to the stars. They 
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proposed to use the Ca n H & K lines to measure the distances to the OB stars. 
They observed that the equivalent widths of the K I and CH lines are less tightly 
related to the parallaxes (distances) of stars. They stated that the K I and CH that 
are present in cloud cores. This observation is in conformance with the McKee & 
Ostriker model because the cold cloud cores have much lower filling factor that the 
cloud corona. 

An analysis of the interstellar Ti II, Ca II and Na I absorption lines was done by 
Welsh et al. (1997) and by Hunter et al. (2006). Column densities of the ionised 
Ti and Ca are well correlated, yielding a conclusion, that both elements reside in 
the same environment. Both Ti n and Ca n interstellar lines analysed by them are 
correlated with distance. 

We want to extend the results presented by Megier et al. (2005) and include 
to the analysis elements observed in the ultraviolet spectral range. Because of the 
difficulties with saturated lines we have used column densities instead of equiva- 
lent widths. The Hubble Space Telescope (HST) spectra were used to calculate 
column densities of Fe II, Ge II, Mg I, Mg II, S I, S II, Si II, Zn II elements. 
Both GHRS (Goddard High Resolution Spectrograph) and STIS (Space Telescope 
Imaging Spectrograph) spectra were used. 

In the next paragraph we present the method used to determine the column den- 
sities. The correlation between column densities and distances is analysed in §3. 
The relationships between interstellar gas-phase abundances and interstellar red- 
dening E(B-V) or f(H2) are also presented in §3. The conclusions are summarised 
in §4. 

2 Column Densities 

The spectra from the ultraviolet spectral range were downloaded from the HST 
Data Archive. The GHRS spectra taken in the FP-SPLIT mode were processed 
with IRAF tasks pojfsets and specalign to achieve the final spectrum. When many 
spectra for one star were available the column density was calculated indepen- 
dently for each spectrum, and an average column density was determined. The 
column densities were derived using the profile fitting technique. The absorption 
lines were fitted by Voigt profiles. The transitions for which the natural dumping 
constant (F) in not known (Mg H 1240 A doublet, Ge H 1237 A, Mg I 1829 A) 
were fitted with a Gauss function. The cloud velocities (v), Doppler broadening 
parameters (b) and column densities (N) for multiple absorption components were 
simultaneously fitted to the observed spectrum. Both lines of magnesium doublet 
(at 2800 A or 1200 A) were also fitted simultaneously - v, b and N were common 
for both Unes in the doublet. The wavelengths, oscillator strengths (f) and natural 
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damping constants (F) were adopted from Morton (2003). 

A convolution with a point spread function (PSF) was also performed. The 
PSF for the GHRS consists of two Gaussian components. The "core" Gaussian 
has a FWHM=1.05 diodes, while the "halo" component has FWHM=5.0 diodes 
(Spitzer & Fitzpatrick, 1993). The relative contribution of the "core" and "halo" 
components into the PSF is wavelength dependent and was interpolated from the 
table given by Cardelli et al. (1990). The Gaussian PSF for the STIS spectrograph 
depends on wavelength, slit and the mode of observations. The tables with FWHM 
for the combination of mode and slit can be found in "STIS Instrument Handbook" 
(Kim Quijano et al., 2003). The FWHM of the Gaussian PSF was wavelength- 
interpolated from these tables. The derived column densities are presented in Table 

m 

The interstellar absorption lines of Fe II, Zn II and S II are often saturated. 
Therefore the column density of these elements was measured only in a few di- 
rections. The results for these three elements are fragmentary, because in the di- 
rections where the absorption lines were saturated, the column densities were not 
determined. 

3 Results 

The column densities of Fe II, Ge II, Mg I, Mg II, S I, S II, Si II, Zn II versus 
parallaxes are plotted on Fig. [T] The parallaxes of stars used in this paper are taken 
from the Hipparcos Catalogue (1997). 

In order to check how the column density (N) changes with distance (d) a 
hyperbola was fitted (see Megier et al. (2005)) 

1 1 

n = — = 

d a- N + b 

in the column density - parallax (tt) plane. We do not use the distance, like Hunter 
et al. (2006) or Welsh et al. (1997), because only for the parallaxes the plus and mi- 
nus errors are equal. Moreover, the parallaxes are sometimes negative. Removing 
such points from the sample leads to the underestimation of distance. The weight 
for each point in the fitting procedure was equal to l/o"J. 

The stars from the Scorpius-Ophiuchus region: HD 141637, HD 143018, HD 
144217, HD 147165, HD 147933, HD 148479 belonging to the Upper Scorpius 
Association (de Zeeuw et al, 1999), and HD 143275, HD 149757 were eliminated 
from the N-;r relation fitting. These are nearby stars lying behind a dense molecular 
cloud. Therefore these directions have very high column densities and deviate from 
the N-;r relation for other stars. The eliminated stars are shown with open circles on 
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Fig. [T] The seven stars from the local interstellar medium (LISM) with n > 20 mas 
were also excluded from the fitting. The results of the column density - parallax 
fitting are following: 

d - (1.23 ± 0.16) • IQ-" ■ N{MgII) + (89 ± 8) [pc] 

d = (1.72 ± 0.25) • lO-^-"* ■ N(SII) + (48 + 10) [pc] 

d = (1.35 ± 0.22) • 10-13 ■ NiSill) + (97 + 11) [pc] 

d = {3.1 ±1.2)- 10"" • mZnll) + (96 ± 38) [pc] 
where column densities (N) are in cm'^, and distances are in parsecs. There are no 
nearby stars with low column density for neutral elements (Mg I, S I) as well as 
for Ge II. Therefore the hyperbola fit gives very high b term, that should be equal 
to the Local Bubble dimension. These fits are shown on Fig. \T\ but these fits are 
probably meaningless. 

Megier et al. (2005) have proposed to use an analogous formula for Ca II to 
determine the distances to OB stars. In our sample of OB stars objects behind 
molecular clouds strongly deviate from the column density - distance relation. 
For example the distance to HD 147933 from the Sco-Oph cloud calculated from 
the Mg II column density is 6 times grater (800 pc vs. 120 pc) than the distance 
determined from the Hipparcos parallax ! Such outliers have f{H2) > 10"^, but 
only some stars with large f{H2) deviate from the N-tt relation. 

The column density of Mg II, Si II and Ge II (dominant ionisation stages) 
normalised to the total hydrogen column density (see Table O versus interstel- 
lar reddening is presented on Fig. |2] The hydrogen-normalised Si II and Mg 
II column density drops for sites with more interstellar dust (grater interstellar 
reddening E(B-V)). For Ge II the drop is not so evident. This can be attributed 
to the lower condensation temperature of germanium (TdGe) = 825 K, versus 
r,(Mg) -r,(50 - 1340 K). 

The physical conditions in the interstellar clouds can be characterised by the 
fractional abundance of molecular hydrogen /(//2) - 2N{H2)/{N{HI) + 2N{H2)). 
The H2 molecule is in equilibrium between formation on dust grains and gas-phase 
destruction in chemical processes or in photodissociation. Some elements , such as 
the noble gas krypton, do not exhibit changes with varying f{H2), because they do 
not condense onto dust grains (Cardelli, 1994). On the other hand elements that are 
easily incorporated into dust grains (Mg, Si) show larger depletions in clouds with 
large f{H2) (molecular clouds) - see Fig. |3] Only a slight decrease of germanium 
abundance was noticed by Cardelli (1994) in a limited sample of directions. In 
our data (Fig. O the increase of germanium depletion is seen in dense, molecular 
clouds. The errors of Ge abundance are large, due to the weakness of the 1337 A 
line. 

The abundances of elements in the dominant ionisation stage are lower than 
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the Solar System abundances for Fe II, Ge II, Mg II and Si II (Fig. gll. The Solar 
System abundances were adopted from Grevesse & Sauval (2000) and are equal to 
{FelH)Q - 3.16 ■ 10"^ 
(Ge/Hh - 3.31 • 10-'^ 
{MglH)Q = 3.80 • 10-5 
{Si/H)Q = 3.59- 10-5 
iS/H)Q - 1.85 • 10-5 
iZn/H)Q - 4.31 • 10"^ 

In the direction to HD 91316 the Zn II abundance is grater than the Solar System 
abundance. Zinc is known to not deplete into dust grains, although the Zn abun- 
dance may be lower than the Solar System abundance {[Zn/S] = log{ZnII/S) - 
log{Zn/S)Q = -0.13, Fitzpatrick & Spitzer (1996)). The S II abundance is grater 
than the Solar System abundance for HD 116658 and HD 68273. Sulphur was 
known to not incorporate into dust grains. Fitzpatrick & Spitzer (1996) and Fitz- 
patrick (1996) have stated that the sulphur abundance is close to the solar abun- 
dance in diffuse clouds. In our data only 2 of 8 directions have the sulphur abun- 
dance close to the solar abundance. The depletion of sulphur reaches [S /H] = 
-1.59 for HD 24398. Such low sulphur abundance leads to doubts if sulphur can 
be used as a tracer of the total hydrogen column density. The low depletion of S 
and Zn is not surprising, since this elements have much lower condensation tem- 
peratures (TciS) = 614 K, TciZn) = 684 K) than other analysed elements (for Fe, 
Mg and Si Tc ~ 1340 K). The sin-like pattern seen on Fig. |4]is caused by changes 
in the interstellar reddening E(B-V) in the analysed directions. 
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star Par Fen Gen Mgl Mgn SI Sn SUI MI 



BD+28 4211 


9.63±1.71 


— 




— 




9.18 + 0.91 


lO" 


5.28 + 1.90 ■ 


10" 


— 




6.93 + 0.28 


lO" 


1.28+0.05 


lO" 


— 




Feige24 


13.44±3.62 


— 




— 




— 




— 




— 




7.53 + 2.64 


10" 


2.76 + 0.15 


10" 


— 




GD246 


— 


— 




— 




— 




— 




— 




— 




9.80 + 0.20 


10" 


— 




HD432 


59.89+0.56 


2.04 + 0.17 


10>2 


— 




— 




3.65 + 0.04 ■ 


10" 


— 




— 




— 




— 




HD 11443 


50.87±0.82 


3.31 +0.77 


10'2 


— 




— 




7.52 + 0.04 ■ 


10" 


— 




— 




— 




— 




HD 18100 


0.88+1.04 


3.39+0.08 


10" 


— 




1.72 + 0.06 


1012 


— 




— 




— 




— 




— 




HD 22049 


310.75±0.85 


— 




— 




— 




3.66 + 0.18 ■ 


10" 


— 




— 




— 




— 




HD 22468 


34.52±0.87 


— 




— 




— 




2.55 + 0.33 ■ 


10" 


— 




— 




— 




— 




HD 23630 


8.87+0.99 


— 




— 




1.07 + 0.03 


10l2 


1.06 + 0.15- 


I0l5 


— 




4.30+1.71 


10" 


8.44+0.11 


10" 


— 




HD 24398 


3.32+0.75 


— 




1.05 + 0.02 


10'2 


— 




2.38 + 0.003 


- 10'5 


— 




7.50 + 2.23 


10" 


— 




— 




HD 24534 


1.21±0.94 


— 




1.42 + 0.11 


10'2 


7.49 ± 0.46 


lO" 


2.71+0.20- 


10l5 


— 




— 




2.01+0.12 


lO" 


— 




HD 24760 


6.06±0.82 


— 




3.65 ± 0.28 


10" 


— 




1.32 + 0.11 - 


10" 


— 




6.14 + 2.27 


10" 


— 




— 




HD 24912 


1.84±0.70 


— 




— 




2.76 ± 0.14 


10" 


3.93 + 0.96 - 


10" 


1.49 + 0.06 


10" 


— 




3.19+0.06 


10" 


— 




HD 28497 


2.07±0.82 


— 




— 




4.86 + 0.11 


1012 


1.11+0.28- 


10" 


— 




— 




1.45+0.05 


10" 


— 




HD 34029 


77.29±0.89 


2.71+0.15 


10l2 


— 




— 




6.77 + 0.14 - 


10" 


— 




— 




— 




— 




HD 35149 


3.39+0.87 


— 




6.83+1.16 


lO'l 


4.36 ± 2.13 


10" 


3.31+0.21 - 


10" 


5.78 ± 0.68 


10" 


1.06 + 0.32 


10" 


2.02 + 0.08 


10" 


7.43+0.09 


10" 


HD 36486 


3.56+0.83 


1.04 + 0.02 


10" 


— 




— 




1.45 + 0.06- 


10" 


— 




— 




1.19 + 0.06 


10" 


— 




HD 36861 


3.09±0.78 


— 




7.41 ± 0.60 


10" 


— 




4.33 + 0.21 - 


10" 


— 




1.94 + 0.71 


10" 


— 




— 




HD 37128 


2.43±0.91 


2.12 + 0.16 


lO" 


2.21 ± 0.54 


lO" 


3.46+1.68 


lO" 


2.14 + 0.08- 


10" 


— 




— 




2.09 + 0.04 


lO" 


— 




HD 38666 


2.52±0.55 


1.90 + 0.02 


10" 


— 




2.91+0.47 


10l2 


1.26 + 0.17- 


10" 


— 




— 




1.49 + 0.05 


lO" 


— 




HD 38771 


4.52±0.77 


— 




3.21 ± 0.75 


10" 


— 




2.20 + 0.14- 


10" 


— 




1.28 + 0.34 


10" 


— 




— 




HD 44743 


6.53±0.66 


— 




— 




3.71+0.22 


10" 


— 




— 




— 




2.40 + 0.20 


10" 


— 




HD 47839 


3.19+0.73 


— 




— 




7.36 + 0.39 


10l2 


2.90 + 0.35 - 


10" 


— 




1.93+0.50 


10" 


2.58 + 0.05 


10" 


6.59 + 0.15 


10" 


HD 57061 


1.02+0.71 


— 




7.22 ± 0.41 


10" 


1.53+0.11 


10" 


5.33 + 0.29- 


10'5 


— 




— 




4.73+0.10 


10" 


1.44 + 0.06 


10" 


HD 61421 


285.93±0.88 


— 




— 




— 




3-65 + 0-33 - 


10" 


— 




— 




— 




— 




HD 62044 


26.68±0.79 


— 




— 




— 




4.78 + 0.20 - 


10" 


— 




— 




— 




— 




HD 62509 


96.74±0.87 


— 




— 




— 




5.33 + 0.65 - 


10" 


— 




— 




— 




— 




HD 68273 


3.88+0.53 


1.06 + 0.01 


10" 










7.81 + 1.49 - 


10" 






1.44 + 0.02 


10" 


1.01 +0.06 


10" 






HD 74455 


1.93+0.57 






5.45+0.51 


10" 


1.08 + 0.03 


10" 


3.07 + 0.01 - 


10" 


2.36 + 0.49 


10" 


3.08 + 0.02 


10" 


3.51 ±0.06 


10" 






HD 89688 


1.86+0.80 


















2.35 ± 0.25 


10" 






3.25 + 0.44 


10" 






HD91316 


0.57+0.82 






2.78+0.73 


10" 


5.22+1.50 


1012 


2.57 + 0.54 - 


10" 










2.57+0.10 


10" 


7.95+0.52 


10" 


HD 100340 


0.70±1.45 






















1.94 + 0.10 


10" 










HD 111812 


10.62±0.90 














3.61+0.03- 


10" 


















HD 116658 


12.44±0.86 














5.51 + 1.73 - 


10" 






2.07 + 0.09 


10" 


1.78 + 1.39 


10" 






HD 119608 


0.77+0.82 










1.27 + 0.24 


10" 


7.97 + 1.03 - 


10" 










4.60 + 0.27 


10" 


1.88+0.15 


10" 


HD 122879 


-0.05+0.76 










1.80 + 0.26 


10" 


1.60 + 0.04- 


10l<5 










1.22+0.05 


lOKi 






HD 141637 


6.25+0.91 






2.45 ± 0.50 


10'2 


1.08 + 0.02 


10" 


7.41+0.19- 


10" 










3.21 +0.18 


10" 


1.38+0.03 


10" 


HD 143018 


7.10+0.84 






6.39 ± 0.94 


10" 


1.39 + 0.06 


10l2 


2.62 + 0.13- 


10" 






2.02 + 0.03 


10" 


2.26 + 0.04 


10" 






HD 143118 


6.61+0.78 














4.47 + 1.45 - 


10" 










4.50 + 0.60 


10" 


1.31+0.12 


10" 


HD 143275 


8.12+0.88 






1.12 + 0.06 


1012 


1.65+0.02 


10" 


5.71+0.25- 


10" 


2.57 + 0.14 


10" 






4,11+0.47 


10" 






HD 144217 


6.15+1.12 










1.05 + 0.65 


10" 






3.20 + 0.47 


10" 






4.16 + 0.09 


10" 


1.50 + 0.11 


10" 


HD 147165 


4.44+0.81 






3.02 + 0.03 


10'2 


1.58 + 0.17 


10" 


8.89 ± 0.24 - 


10" 


3.95+0.46 


10" 






4.93 + 0.29 


10" 






HD 147933 


8.27+1.18 






2.20 + 0.46 


1012 


3.12 + 0.26 


10" 


5.79 + 0.64- 


10'5 


3.52 + 0.63 


10" 






3.60 ± 0.20 


10" 






HD 148479 


5.40+1.68 






1.92 + 0.26 


1012 


4.32 + 0.63 


10" 


1.67 + 0.08 - 


lOlo 


2.89+1.46 


10" 






6.37 + 0.06 


10" 






HD 149499 


26.94+1.88 






















1.13 + 0.17 


10" 


7.39 + 0.44 


10" 









star 


Par 


Fen 


Gen 




Mgl 




Mgn 




SI 




sn 


SUI 




Znn 




HD 149757 


7.12±0J1 


— 


8.98 + 2.38 


lO" 


2.38 + 2.11 


lo'J 


2.61 + 0.22 


10" 


3.42 + 0.46 


lO" 


— 


2.10 + 0.04 


lO'* 


— 




HD 154368 


273+0.96 




3.68 + 1.58 


lO'^ 






1.04 + 0.11 


lo" 


3.41 + 0.92 


lO" 




5.81 ±0.61 


lo" 






HD 158926 


4.64+0.90 








1.05 + 0.29 


10" 












1.56 + 041 


10" 






HD 160578 


7.03+0.73 








8.96 + 2.33 


lO" 


4.39 + 0.57 


10" 








9.19 + 0.16 


10" 


4.21+0.26 


1012 


HD 195965 


1.91+0.58 




1.39 + 0.28 


10'2 


4.67+0.84 


10" 


7.96 ± 0.80 


10" 


1.73 + 0.17 


10" 




5.56+0.14 


10" 






HD 198478 


1.45+0.55 




3.19 + 0.25 


10'2 






8.17 + 0.80 


10" 
















HD 201345 


0.61+0.77 












9.14 + 0.80 


10" 
















HD 202904 


3.62+0.56 








2.90 + 0.33 


1012 


1.58 + 0.24 


10" 








1.21+0.04 


10" 


6.09 + 0.15 


1012 


HD 203374 


-0.70+0.70 




2.76 + 0.31 


10'2 


1.11+0.06 


lO" 


1.22 + 0.05 


10'« 


4.45 + 0.11 


10" 




1.41 +0.19 


1016 






HD 203664 


2.23+1.08 








1.26 + 0.67 


10'^ 


2.93 ± 0.59 


10'5 








4.37 ± 0.23 


10" 






HD 206267 


2.78+079 




2.92 + 0.21 


10'2 


1.80 + 0.09 


lO" 


1.05 + 0.05 


10"> 








1.26 + 0.15 


10'<' 






HD 207198 


1.62+0.48 












1.15 + 0.05 


ioi<> 
















HD 207538 


0.30+0.62 












1.06 + 0.06 


lOKS 
















HD 209339 


0.14+0.57 




2.31+0.17 


1012 


1.02 + 0.11 


10" 


1.06 + 0.02 


lOlf' 


2.56 + 0.13 


10" 




9.43 + 1.43 


10" 






HD 210839 


1.98+0.46 




3.12 + 0.92 


1012 


1.43 + 0.06 


10" 


1.13 + 0.06 


10^<> 


4.34 + 0.15 


10" 




1.05+0.07 








HD 212571 


2.96+0.71 








1.29 + 0.77 


10" 


2.50 + 0.30 


10" 


2.88 + 0.69 


1012 




2.23+0.04 


10" 






HD 219188 


0.94+1.16 








1.36 + 0.18 


10" 


3.95 ± 0.20 


10'5 








3.54 + 0.18 


10" 






HZ 43 


31.26+8.33 












1.84 + 0.26 


10'2 

















Table 1: Parallaxes [mas] and column densities [cm ^] for the analysed directions. 



4 Conclusions 



• The column density of ionised elements Mg II, S II, Si II, Zn II correlates 
with the distance. 

• Estimation of a distance from the column density can lead to large overesti- 
mation of the distance for star lying behind molecular clouds. 

• Sulphur is depleted in comparison to the Solar System abundance. Sulphur 
should not be used as a tracer of the hydrogen column density. 

• Germanium abundance drops with increasing fractional abundance of molec- 
ular hydrogen. 
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Table 2: Interstellar reddening and hydrogen column density for target stars. Ref- 
erences: 1 - Bohlin et al. (1978) ; 2 - Diplas & Savage (1994) ; 3 - Jenkins et al. 
(1986) ; 4 - Lacour et al (2005) ; 5 - Rachford et al (2002). 
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Figure 1: Column densities [cm ^] versus parallaxes [mas]. The grey open circles 
represents stars from the Scorpius - Ophiuchus region. 
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Figure 2: Column densities of Mg 11, Si 11 and Ge 11 normalised to the total hydro- 
gen column density (logarithmic scale) versus the interstellar reddening E(B-V). 
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Figure 3: Column densities of Mg 11, Si 11 and Ge 11 normalised to the total hy- 
drogen column density versus fractional abundance of molecular hydrogen f(H2). 
Both axis use logarithmic scale. 
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Figure 4: Black dots show the abundance of the dominant ionization stage rela- 
tive to the Solar System abundance [dex]. The grey lines represent the interstellar 
reddening E(B-V). The stars are ordered by HD numbers. 



14 



